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We study current-induced dynamics of a magnetic domain wall by solving a time-dependent 
Schrodinger equation combined with Landau-Lifshitz-Gilbert equation in a one-dimensional 
electron system coupled to localized spins. Two types of domain-wall motions are observed 
depending on the hard-axis anisotropy, K±, of the localized spin system. For small values of 
K±, the magnetic domain wall shows a streaming motion driven by spin transfer. In contrast, 
for large values of K±, a stick-slip motion driven by momentum transfer is obtained. We clarify 
the origin of these characters of domain- wall motions in terms of the dynamics of one-particle 
energy levels and distribution functions. 
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A magnetic domain wall (DW) is a twisted structure 
of localized spins, which separates the magnetic domains 
with different polarization in a ferromagnet. Dynamics of 
topological defects, such as DWs and vortices, has been 
one of the central issues in condensed matter physics. The 
lowest excitations of polyacetylene are described as soli- 
tons, which determine the conductive properties of the 
system. 1 The critical current of type-II superconductors 
is sensitive to the pinning of vortices. 2 Recently, current- 
induced dynamics of magnetic DWs attracts consider- 
able attention for its importance as a non-equilibrium 
many-body problem and its potential for technical appli- 
cation. 3 5 In particular, promising proposals have been 
made for the application to non- volatile memory. 6 Nano- 
scale fabrication of the memory device may be possible, 
by controlling the DW motion by an electric current. For 
simplicity, let us consider a situation where there is a sin- 
gle DW separating the -f z-spin and the — z-spin regions 
(Fig. 1), which is coupled to electrons ferromagnetically. 
Then suppose that an electron wave packet with -f z-spin 
is injected from the left-hand side. In this case, one can 
imagine two mechanisms for the dynamics of the DW. 7 ' 8 
(i)The spin of electron wave packet staying in the -f z-spin 
region is aligned in the +z direction due to the ferro- 
magnetic coupling to the localized spins. If this electron 
is carried adiabatically to the — z-spin region, the direc- 
tion of electron spin gradually changes and ends up in 
—z. Due to the conservation of angular momentum, the 
total z-axis spin of localized spins increases by fi, which 
leads to the enlargement (shrinking) of the +z (— z)-spin 
region, resulting in a shift of the DW. This adiabatic pro- 
cess is called "spin transfer." (ii)On the other hand, if the 
electron is transferred to the — z-spin region with a finite 
velocity, the electron is reflected by the DW with a finite 
probability. Through this reflection, the DW acquires a 
finite momentum as a result of recoil effect. This pro- 
cess is called "momentum transfer," and originates from 
the non-adiabatic motion of the conduction electrons, in 



contrast to the "spin transfer." 

The character of DW motion becomes qualitatively dif- 
ferent depending on which mechanism is dominant. So 
far, extensive theoretical studies 8 14 have shown that the 
"spin transfer" is the dominant mechanism for the DW 
motion in the experimental situations. Actually, many 
experimental results have been explained in this mech- 
anism. 15 19 On the other hand, the relevance of "mo- 
mentum transfer" has been also reported in some exper- 
iments. 20,21 The dynamics of thin DW also draws atten- 
tion, 22 in which more careful treatment is needed for the 
effect of non-adiabaticity of conduction electrons, i.e., the 
effect of "momentum transfer." 

Most of the theoretical studies have implicitly assumed 
that the conduction electrons obey a static distribution 
function slightly off the equilibrium state. Under such as- 
sumption, the non-adiabaticity of conduction electrons 
is not correctly taken into account. In particular, it is 
expected that the effect of "momentum transfer" is un- 
derestimated. Therefore, in this paper, we will treat the 
dynamics of conduction electrons and a DW on an equal 
footing, and discuss how the DW motion is affected by 
the strong non-adiabaticity of the conduction electron 
system. 

Our model consists of a one-dimensional electron sys- 
tem under an electric field, which is coupled to a classical 
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Fig. 1. (Color online) The schematic figure of a DW formed by 
localized spins: X and </>o are the displacement and the rotating 
angle of DW, respectively. 
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localized spin system, 

n(t) = -t c J2[e im 4 +1(T c ia + h.c.} 

~ ^sd ^ ] Si(t) ' cj a (7 'aa'Cicr' ' . 



(i) 



Here, ^(cj^) is an annihilation (creation) operator of the 
conduction electron at the zth site and spin a under the 
periodic boundary condition, Ci+Na — c ia- Si(t) denotes 
classical localized spins at the sites z, which constitute 
DWs with an easy z-axis and a hard y-axis. The electric 
field is applied via a time-dependent AB flux, &(t). In 
order to apply a direct electric field, we set $>(t) = 
where V is the voltage drop in the circuit. 

We study the dynamics of state vector of conduction 
electrons, and the localized spins, Si(t), simul- 

taneously, is determined by the time-dependent 

Schrodinger equation, \hf t \^{t)) = H(t) At the 

same time, we determine Si(t) by solving the semiclas- 
sical Landau-Lifshitz-Gilbert (LLG) equations, 



Si = H em xSi- a(Si/S) x Si 



(2) 



The first term of eq. (2) represents the torque which leads 
to the precessional motion of the localized spin around 
the effective magnetic field, iieffi- The second term cor- 
responds to the Gilbert damping, which is essential to 
the relaxation of Si parallel to —H^. Here, 



H e ffi = —J(Si-i + Si+i) - KS iz e z 
+ Kj_Si y e y — J s d (<Ti) , 



(3) 



with the exchange coupling between localized 
spins J, the easy-axis anisotropy the hard- 

axis anisotropy K±, and the electron spin (<Ti(t)) 

= 4*°™'**' at the ith site - 

In order to solve the time-dependent Schrodinger 

equation numerically, we use the Crank-Nicholson's 
method, 23 \9(t + At)) = exp[-^ fi +At H(t')dt'] 

- iT§7§W§T^7§^ which ensures the norm 

conservation of the state vector. For the LLG equations, 
we use the implicit fourth-order Runge-Kutta method. 

We set the initial configuration (t = 0) 
of the localized spins to be the static solu- 
tion of the LLG equations, eq. (2), with two 
DWs: 5,(0) = 5{l/cosh[^f^],0,tanh[^if^]} 
for the sites 1 < i < N/2, and S;(0) = 
[-Si.jv^xWA-W^'-^-WO)] for N/2 + 1 < 
i < N. Here, i c i(0) = (N + 2)/4 is the initial coordinate 
of the first DW center, and W is the width of the 
DWs. Note that an even number of DWs are required 
for smooth connection at the boundary because of 
the periodic boundary condition. Since the relation 
Si+N/2(t) = [Six(i),Siy(i), -S iz (t)} holds within the 
simulation time, we only focus on the dynamics of the 
first DW. We assume the initial state vector, |^(0)), to 
be the ground state of Hamiltonian, eq. (1), with the 
initial configuration, Si(Q). 

In order to characterize the DW motion, we define the 
displacement of DW center, X(t), and the rotating an- 
gle of the localized spins in the DW center around z- 
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Fig. 2. (Color online) The time dependence of (a)the displace- 
ment X, and (b)the rotating angle 4>q of the DW for several 
values of the hard-axis anisotropy K±. The horizontal line of (b) 
corresponds to 4>q = 7r/4. The scale of the longitudinal axis is 
expanded in the insets of (a) and (b). 
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Fig. 3. (Color online) The time dependence of the velocity, X, 
and the force exerted by the electrons, F, for K± = 1. The 
amplitude of F is plotted in arbitrary unit. 



axis, 0o (t) , as shown in Fig. 1. We will calculate the z- 
component of spin torque, r z (t) = J s d SieistDwI^W x 
{(Ti(t))) z . Since the first term in the right-hand side of eq. 
(2) involves the term proportional to r z (t), Si Z increases 
with r z (t). This means that r z (t) transfers z-axis angular 
momentum from conduction electrons to localized spins, 
and is essential to the "spin transfer" mechanism. We 
also calculate the force, F(t) = — J s d SieistDW W ' 
(<7i(t)), which is the time variation of the momentum 
passed from the electron system to the localized spin 
system. Note that only if the DW has mass, F(t) can 
be identified as a 'real' force exerted by the conduction 
electrons on the DW and is relevant to the "momentum 
transfer" mechanism. As shown below and derived in Ref. 
[8], both r z (t) and F(t) determine the dynamics of DW. 
In the following, we assume that t c and h/t c are the units 
of energy and time, respectively, and that h = e = S = 1. 

Here, we will show our results. The time evolution of 
X(t) and (j>o(i) are shown in Figs. 2(a) and 2(b), respec- 
tively, for several values of K±. We set N = 60, W = 5, 
J = 2.59, K = 0.1, a = 0.02, J sd = 0.3, V = it x 10~ 4 , 
and the electron density n = 12 /N = 0.2. As shown 
in Figs. 2(a) and 2(b), the DW motion qualitatively 
changes, depending on the value of K±. 
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Firstly, for small values of K± (K± = 0.01, 0.02, 
0.025), X(t) and (/>o(t) increase smoothly with t as shown 
in the insets of Fig. 2. The constant increase of (j)o(t) is 
caused by the spin torque r z (t) with the help of Gilbert 
damping. 24 This streaming motion of the DW accom- 
panied by the continuous increase of 4>o(t) means domi- 
nance of the "spin transfer" mechanism. 8 

Secondly, for large values of K± (K± = 0.03, 0.05, 0.1, 
1), step-like increase of X(t) with t is observed in Fig. 
2(a). In this case, the rotating angle <po(t) is saturated 
at a finite value in a longer time scale in contrast to the 
case of smaller K±. We plot the velocity X(t) and the 
force F(t) for K± = 1 in Fig. 3. Remarkably, X{t) and 
F(t) show peaks at the same time. This perfect corre- 
spondence clearly shows that the DW motion is driven 
by the force given by conduction electrons, i.e., the "mo- 
mentum transfer" is the dominant driving mechanism in 
this region. We will later discuss the microscopic origin 
of the periodic peak structure of F(t). 

Here, we discuss the origin of crossover in driving 
mechanism between "spin transfer" and "momentum 
transfer." The crossover results from the competition be- 
tween the spin torque r z (t) and the hard-axis anisotropy 
K±. K± tends to fix (/>o(t) ~ 0, while r z (t) increases <po(t) 
as we mentioned in the footnote [24]. The origin of the 
crossover behavior can be interpreted using the effective 
equations of motion for X(t) and 0o(t), 8 

2(-fa + aX/W) =F, (4) 
2{X + aW<j> ) = r z - K± W sin 2^ , (5) 

where irrelevant coefficients are omitted. For r z (t) > 
K±W, since the right-hand side of eq. (5) is always fi- 
nite, constant supply of the spin torque r z (t) ^> aWF(t) 
keeps both X(t) and (/>o(t) ~ aX(t)/W increasing. This 
is why the DW motion for smaller K± is driven by the 
"spin transfer" (the insets in Fig. 2). In contrast, for 
r z (t) < K±W, the effect of the spin torque is canceled by 
the contribution of large K±, and then 4>o(t) is saturated 
at the angle (j) s = - sin -1 ( K T ^ W ) (< 7r/4) in a longer time 
scale (K ± > 0.03 in Fig. 2(b)). This is the so-called in- 
trinsic pinning effect for the "spin transfer" mechanism. 
In the case of r z (t) <C K±W, since 0o(£) <C 7r/4 and 

mass of the DW is well-defined as Mow = ^k^v? > ^ ne 
"momentum transfer" from the conduction electrons be- 
comes relevant to the DW motion. 

With r z (t) <C K±W in mind, substituting (j>o(t) = 
(j) s + 5(/)(t) into eqs. (4) and (5), we obtain the equations, 
2(— 50 + olX/W) = F, and X + aW5<j) = -K±W5</>. 
The solutions are estimated as 5(j)(t) = — 2(i+a 2 ) F(t) 

and X(t) = 2a[i+Q2) ^ft) ? including a numerical factor, 
< c < l. 25 In our calculation for K± = 1, c ~ 0.95 is 
obtained by numerically averaging the ratio of X(t) and 
F(t) over the interval: 4 x 10 4 < t < 1.2 x 10 5 . Since 
a is small in general, the contribution of the force to 
the DW motion is remarkably enhanced, being inversely 
proportional to a. 26 The fluctuation of 4>o(t) around the 
saturation angle </> s is crucial for the enhancement of the 
DW motion: X oc a -1 . 

Our calculation shows the absence of the perfect in- 



trinsic pinning in a large K±_ region. While (j>o(t) is sat- 
urated for K± > 0.03, the DW still shows a step-like 
motion. Absence of the intrinsic pinning has also been 
reported in the numerical study of ref. [12]. However, 
they attribute the origin of DW motion to the non-local 
deformation of localized spins, which is different from our 
results. In our calculation, deformation of localized spins 
is negligible. 
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Fig. 4. (Color online) One-particle energy levels {£^(0)} for var- 
ious J S( j values. The energy levels split into two bands. For 
<^sd ^ fi/2, the two bands overlap at |^(0)| < —J s d + ^/2 
where spin and momentum degrees of freedom are entangled, 
and a complicated energy level structure is formed. 
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Fig. 5. (Color online) The time dependence of the one-particle 
energy levels {Ei(t)} (thin lines) and the one-particle distribu- 
tion functions {| (^(t)|^(t)) | 2 } (color plots) of the electron sys- 
tem for K± = (a) 1 and (b) 0.01. Note that the time scale in 
(b) is much larger than that in (a). The color of the thick lines 
represents the occupation probability of each level. The dot and 
the arrow in (a) denote an elastic backscattering and an inelastic 
scattering, 29 respectively. 



In order to understand the periodic step-like motion 
of the DW from a microscopic point of view, we will 
study the dynamics of electron state vector, \Sfr(t)). Be- 
fore studying the dynamics, let us look at the one-particle 
energy levels at t = 0, i.e., when there is a stationary 
DW. Figure 4 shows the energy levels {£^(0)} of the 
Hamiltonian H(t) (eq. (1)) at t = for various values 
of J s d- At J s d = 0, the energy levels are simply those of 
one-dimensional tight-binding chain with a band width, 
Q = 4. Each level is at least four- fold degenerate in terms 
of the reversal of momentum and spin. With increasing 
J s d, the spin degeneracy is lifted, and the energy lev- 
els gradually split into two bands. The lower (higher) 
band, -(+) J sd - < £;(0) < -(+) J sd + fi/2, is com- 
posed of the states where the spins of conduction elec- 
trons are aligned parallel (antiparallel) to the localized 
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spins. For not too large J s d (J s d ^ ^/2), the two bands 
have an overlapping region, |£^(0)| < — J s &+£l/2. Within 
each band, levels are labeled with different momentum 
except for the overlapping region, where neither spin nor 
momentum serve as a good quantum number, i.e., the 
eigenstates are composed of entangled spin and momen- 
tum degrees of freedom. This "entangled" energy range 
will play an important role in the dynamics in the larger 
K± region, as we will show below. 

To characterize the dynamics of | we diagonal- 

ize the Hamiltonian TL(t) (eq. (1)) at each t, and obtain 
the eigenvalues {Ei(t)} and eigenvectors {!&(£))}• We 
then plot the time evolution of one-particle energy lev- 
els {Ei(t)} and distribution functions {| | 2 } 
in Fig. 5, corresponding to the DW motions for K± = 1 
and K± = 0.01. First, let us discuss the case of large K± 
(Fig. 5 (a)). In general, the time dependence of {Ei(t)} 
comes from the AB flux, and the localized spins, 

Si(t) in eq. (1). For K± = 1, the dynamics of {Ei(t)} is 
mostly determined by since Si(t) is constant except 
for the moments at which X(t) shows peaks (Fig. 3). §(t) 
accelerates the electron momentum with an acceleration 
rate: dk/dt ~ V/N. As to the distribution functions, only 
the states below the Fermi energy (epo ~ —1.98) are oc- 
cupied at t = 0. Therefore, for small t, all the occupied 
states are within the range Ei(t) < J S( ± — fl/2. In this en- 
ergy range, each state continues to be accelerated with- 
out scattering, since momentum serves as a good quan- 
tum number. 27 On the other hand, when Ei(t) enter the 
"entangled" energy range, \Ei(t)\ < — J s d + ^/2, an elas- 
tic backscattering with spin flip occurs, and as a result, 
the DW acquires a recoil momentum. 28 The backscatter- 
ing causes rapid oscillation of the localized spins (see the 
oscillating behavior of X(t) just after the peaks in Fig. 
3), which leads to the dissipation (arrow in Fig. 5(a)). 29 

With the difference of momentum between the suc- 
cessive energy levels at t — 0, Sk ~ 2tt/N, the in- 
terval between the backscatterings St is estimated as 
St ~ Sk/ (dk/dt) ~ 2tt/V. Using V = tt x KT 4 , we have 
St = 2 x 10 4 , consistent with Fig. 5(a), where backscat- 
terings happen at time t e b ~ 5 x 10 4 , 7 x 10 4 , 9 x 10 4 
Remarkably, the backscatterings occur exactly at the mo- 
ments at which F(t) has peaks (Fig. 3). This correspon- 
dence shows that the backscatterings in the "entangled" 
energy range is the origin of the periodic step-like DW 
motion. Actually, the most part of DW shift Z\X ave ~ 65 
(Fig. 2(a)), which is a time-averaged value of X be- 
tween t e b — St/2 < t < £ e b + St/2, can be explained 
by the "momentum transfer" at the elastic backscatter- 

ing, AX eh = Z[ZS] Ak ~ 62 > where Ak ~ mk l 2 
is a momentum shift per DW at t^ estimated in Fig. 
5(a). It is tempting to associate the step-like DW mo- 
tion with the stick-slip motion: 30 "stick" caused by the 
pinning due to the hard-axis anisotropy K±, and "slip" 
caused by backscatterings of the electrons. The stick-slip 
motion vanishes, if the lifetime of conduction electrons 
is shorter than 2tt/V. However, it may be observed in a 
clean mesoscopic sample. 

Next, we discuss the case of small K±, i.e., K± = 0.01 
(Fig. 5(b)). The qualitative feature of level dynamics is 



distinct from that for K± = 1. As soon as the electric 
field is applied at t = 0, the energy levels of the lower 
band start to reconstruct themselves so that each level 
has almost time-independent energy, except for small os- 
cillations. As a result, the system falls into a nonequilib- 
rium steady state, where the inter-level transition rarely 
happens. The distribution function of the system can be 
well approximated by a shifted Fermi distribution func- 
tion, with the total momentum proportional to the ap- 
plied electric field. 8 ' 13 

In summary, we have studied the current-induced 
dynamics of a DW by solving the time-dependent 
Schrodinger equation and the LLG equations simultane- 
ously for a one-dimensional tight-binding chain coupled 
to a localized spin system. We observe two types of DW 
motions depending on the hard-axis anisotropy of the lo- 
calized spin system, K±. For smaller K±, the DW shows 
a streaming motion driven by the spin transfer. In con- 
trast, for larger K±, we obtain a stick-slip motion driven 
by the momentum transfer. To our knowledge, our work 
provides the first numerical realization of the crossover 
between the two regimes. We have also interpreted the 
qualitative features of the DW motion in both regimes, 
in terms of the dynamics of one-particle energy levels and 
distribution functions. In particular, we found that the 
stick-slip motion at larger K± comes from the backscat- 
terings in the "entangled" energy range. Our results in- 
dicate that the DW motion at larger K±_ is influenced by 
the microscopic energy level structure, which has been 
neglected in most of the previous works. The competition 
between the stick-slip motion and an extrinsic pinning or 
impurities is an interesting future problem. 19 22 
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